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^\. • Abstract. The recently discovered Submillimeter Extragalactic Background (submm 

On ' EB) has opened new perspectives on our understanding of galaxy evolution. We 

^^ , detail in this paper one of the major cosmological consequences of this discovery, 

the submrn luminosity density history of the Universe. 
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1. The Extragalactic Background (EB) 

^ . Our knowledge of the early epochs of galaxies has recently increased 

OO I thanks to the observational evidence provided by UV/Vis/Near-IR, far- 

^*^ • IR and submm surveys of high-redshift objects. In a consistent picture, 

^-^ . the epoch of galaxy formation can also be observed by its imprint 

on the background radiation which is produced by the line-of-sight 
^^ • accumulation of all extragalactic sources. 

The search for the UV, optical, near-IR and mid-IR EB (from 2000 A 



i-Q ' to 15 /xm) obtained by summing up the contributions of galaxies using 



Dh' 



number counts currently gives only lower limits. However, in the optical 



O I and near-IR the flatening of the faint counts obtained in the Hubble 

-4— » ' Deep Field suggests that we are now close to convergence. Thus we 

C^ , already have a good determination of the background in the optical 

and near-IR. 

At longer wavelengths (A >100 /im), a submm EB was predicted to 
r> I exist as early as the 1970s. It has been detected for the first time 

C^ ■ only in 1996 by Puget et al. thanks to COBE data. Its detection, 

now confirmed by several other studies (Fixsen et al. 1998; Hauser 
et al. 1998; Lagache et al. 1999) has provided new perspectives for 
understanding galaxy evolution. The amount of energy in the submm 
EB is rather large. It exceeds the predictions based on extrapolations 
of the starburst galaxies as seen in the IRAS deep surveys by about a 
factor of 3 (Franceschini et al. 1994). Thus, it requires other sources, for 
example, spheroidal systems radiating mostly in the far-IR during their 
early evolution, or a high formation rate of massive stars at the early 
stages in the evolution of elliptical galaxies. Whatever the nature of the 
emitting sources, the background implies a very strong cosmological 
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evolution from the local galaxies to the more distant ones (Guiderdoni 
et al. 1998; Franceschini et al. 1998), with an enhancement of the IR 
emission at early stages of the galaxy evolution. This very strong evolu- 
tion is also supported by the comparaison between the energy contained 
in the optical and submm EB. The sumtotal of the fluxes of the EB 
from 6 to 1000 fim is about 3.3 10~® W m~^ sr~^ compared to 2 10~^ 
W m~^ sr~^ between the UV band and 6 /im. For the local galaxies, 
the energy ratio between IR and optical wavelengths is much smaller 
than that measured from the background (about 0.4 compared to 1.6). 
This simply reflects a strong change of galaxy properties between the 
local Universe and the more distant Universe. This change must affect 
the Star Formation (SF) history derived only from optical and near- 
IR observations of high-redshift objects. A significant fraction of SF 
might be hidden in heavily extinguished galaxies. These galaxies can 
be totally missed by optical surveys. 

2. Star formation history 

The history of the cosmic Star Formation Rate (SFR) can be derived 
from deep optical surveys assuming that (1) the stellar Initial Mass 
function is universal, (2) the far-UV light is proportionnal to the SFR 
and (3) extinction is negligible. The presence of dust which absorbes 
the UV starlight makes this last assumption not valid. The corrections 
needed to account for the extinction are rather uncertain and there 
is much controversy about the value of this correction. Moreover, the 
SFR deduced from optical surveys can be underestimated if there is a 
significant population of objects so obscured that they are not detected 
in these surveys. A way to avoid this dust contamination is to derive 
the SFR from the IR/submm surveys. However, so far, the catalogues 
of faint submm sources with reliable redshifts are not large enough 
to reconstruct the history of the SFR. But strong constraints can be 
provided by the inversion of the submm EB. 

2.1. Inversion of the submm EB 

For a Universe following the Robertson- Walker metric, the intensity of 
the submm EB (1^) can be simply related to the number of sources per 
Mpc'^ (Nz) having a luminosity L^g using the expression: 

vly = --\ / NzLy^^—dvQ b{vQ--— — )dz (1) 

4vr 7^=0 Jvo=u dz 1 + z 

Choosing a spectral energy distribution Lj^,, and a cosmological model 
(i.e dt/dz), the inversion of Eq. 1 will give A^^. The variation of the 
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luminosity density as a function of redshift ^{z), is then derived as 
follows: 

>f{z) =N, I L^^duo Lq/Mpc^ (2) 

Using such a formalism, the production rate of the far- infrared radiation 
as a function of redshift can be deduced using only one assumption: the 
spectrum of the galaxy L^^,. Since the submm EB seems to be domi- 
nated by luminous IR galaxies, as indicated by the ISO (Aussel, 1998) 
and SCUBA (e.g. Lilly et al. 1999) results, we choose a spectral energy 
distribution L^^ typical of starburst galaxies (Maffei, 1994; Guiderdoni 
et al. 1998). To determine cp{z) (Eq. 2), we choose to explore the cases 
given by the combinations of: (1) three infrared galaxy luminosities: 3 
10^^ L0, 5 10^^ Lq and 3 10^'^ L© (2) two values for the dust spectral 
index (-1.7 and -2) and (3) three cosmological models defined by the 
set of parameters h, Dq and Q\ (h=0.65, r2o=0.3, i^A=0.7; h=0.65, 
Oo=0.3, rJA=0, and h=0.65 ^^0=1, ^A=0) which fixes the dt/dz. The 
basic algorithm for finding N^ is based on Monte Carlo simulations. N^ 
is sampled at a few redshift values and linearly interpolated between 
these values for computing the term on the right of Eq. 1. It has been 
assumed that beyond z=13, N2=0. The best evaluation (minimum x^) 
has been obtained by exploring a wide range of randomly distributed 
values of N^ . Error bars are estimated by keeping the computed submm 
EB (1) within the submm EB uncertainties at each sampled frequency 
and (2) greater than the lower submm EB limit at 850 ^m from Barger 
et al. (1999). By using an iterative method for progressively reducing 
the range of explored values, we reach convergence for each case studied 
with a reasonable number of hits (~ 50000). For the submm EB {I^, in 
Eq. 1), we use the FIRAS and the DIRBE 140 and 240 ^m determi- 
nations of Lagache et al. (1999). (p{z) at z=0 has been fixed following 
Soifer & Neugeubauer (1991): for h=0.65, we have ^p{z = 0)=1.07 10^ 
Le/Mpc^. 

2.2. Main result 

Results for (p{z) are presented in Fig. la for a fixed cosmological model 
and three different luminosities. We note that the combination of all 
parameter sets gives results consistent within a factor of two. As was 
expected, constraints on (p{z) are very weak below redshift 1 (no EB 
values below 140 nm) and above 4 (very low signal to noise ratio of the 
FIRAS submm EB above 800 //m). Between redshifts 1.5 and 3.5, the 
submm EB gives strong constraints. The luminosity density is about 
10 times higher at z=1.4 than at z=0 and it is nearly constant up to 
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Figure 1, (a) Luminosity density distribution as derived from the submm EB for a 
given cosmological model (h=0.65, f2o=0.3, r2A=0.7) and dust spectral index (q = 2) 
with three different luminosities (stars: 3 10^^ L©, diamonds: 5 10^^ Lq, triangles: 
3 10^" Lq). The error bars represent the uncertainties derived for the star points. 
Also represented is the analytical function that reproduces the case L=5 10^^ Lq, 
^0=1 and Oa^O. The analytical function is io°^'+6^'+"+d with a=-2.42 10"^ 
b=5.78 10"^, c=1.36 and d=8.03. (b) SFR derived from the submm EB (analytical 
function, continuous line) and UV/Vis/Near-IR observations (diamonds: Lilly et al. 
1996; triangles: ConoUy et al. 1997; squares: Madau et al. 1996 and crosses: Steidel 
et al. 1999) 

redshift 4. Data points derived from optical observations (Fig. Ib)^ are 
lower by a factor of about 5 than our determination. This clearly shows 
that the optically-derived SFR is severely underestimated due to the 
presence of dust. 
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To compare the star formation rate with the luminosity density, we use: 
fy^ (Guiderdoni et al. 1998) 



'iQyr 
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